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NASA/BLM APT FINAL REPORT 
VOLUME It EXECUTIVE SUMMARY 


1.0 INTRODUCTION . 

1.1 Organization of Final Report . 

The final report for the Phase II activities on the 
NASA/BLM Applications Pilot Test (APT) is divided into three 
stand-alone documents. Volume I is the Executive Summary, and 
as such describes the Phase II project in general terms on 
objectives, procedures and results. Volume ll. Technology 
Demonstration, describes the methodology and results of the 
vegetation mapping and inventory activities for the Arizona Test 
Site. Volume III, Technology Transfer, details the training 
courses presented to BLM personnel on the Integration and use 
of remotely sensed data for mapping and inventory requirements 
as implemented during the Phase II project. 

1.2 Background and Direction of APT . 

The basis of the current APT program jointly being 
carried out by NASA and BLM was a remote sensing research 
project contracted by BLM to the University of California, 
Berkeley (Colwell 1975). That project, popularly referred to 
as the Susanville Program, took place in BLM's Susanvllle 
District in Northeastern California and Northwestern Nevada in 
the 1973-74 time period. Although the program was concerned 
with remote sensing in its most general terms, a number of 
specific applications of Landsat and aerial photography at 
various scales -.were Investigated. An analysis of output 
products from that study concluded that certain information 
extracted from satellite data and supported with other stages 
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1 . 2 --Continued . 


of data could be used to prepare unit resource analyses and 
thence for management framework and activity planning. 
Specifically, the Susanville Program developed vegetation 
cover maps and production estimates. 

As a result of this experience personnel from BLM 
and NASA met in May 1975 to explore the basis for a joint 
remote sensing project. At BLM's suggestion, an inventory of 
wildland vegetation was selected since vegetation plays a role 
in all BLM resource activities i^ands, wildlife, minerals, 
recreation, watershed, range, ^.id forestry) . By February 1976 a 
project had been prepared and in November 1976 a memorandum of 
understanding between BLM and NASA was signed. The NASA/BLM 
program was initiated as an Application System Verification and 
Transfer (ASVT) and was later change^ to an Application Project 
Test (APT) . The program undertaken was a phased three year 
project which encompassed three test sites: in Alaska, in 
Arizona, and in Idaho respectively. The program was designed to 
initially expose the user agency (BLM) to technological 
concepts and procedures in remote sensing and thereafter to 
have BLM assume an ever increasing responsibility for the 
conduct of the program. At its conclusion total responsibility 
for planning and implementation will reside with BLM. The 
first (Alaska) phase was carried out by NASA working through a 
contractor (ESL Incorporated) and was for the purpose of 
vegetation cover mapping in a test site representative of a 
northern spruce tundraCibiome. ] During this phase BLM assisted 
with the field work, witlT-the' digital processing and analysis 
and with monitoring the contract. The second (Arizona) phase 
of the APT program has included both vegetation cover mapping 
and range, woodland, and forest production estimation. The 
former task (cover mapping) is being carried out by BLM. The 
latter task, to which this report is primarily directed, is being 
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carried out by ESL Incorporated under contract with NASA. 

The Phase II site is located in extreme Northwestern Arisona 
and is typical of the southwest desert and upland communities. 

The third (Idaho) phase is to include both vegetation cover 
mapping and range forage production estimation, both tasks in 
this case being carried out by BLM. The Phase III site lies 
in Southwestern Idaho and is typical of a Great Basin sagebrush/ 
grassland community. 

A joint planning session including NASA/ JSC, BLM/DSC 
and Arizona Strip District, and ESL and their subcontractors 
was conducted on October 3-4, 1978 prior to beginning the 
contract effort for Phase II in Arizona. The purpose of the 
session was to finalize project philosophy, objectives and 
schedules and workshop (training) content, attendees, location 
and schedule. A similar session at the beginning of Phase I in 
Alaska was very helpful as a beginning point in that the 
participants were able to exchange views and have an overall 
view of the tasks to be accomplished. 

ESL has prepared a flow chart with a schedule of 
tasks and tentative dates. Using this as a starting point, 
a final flow chart of events for the entire Phase II effort 
was developed. A copy is included as Appendix I. Although 
changes were required, the flow chart was an invaluable tool 
as a reference point for identifying the impact of changes on 
other events and schedules. 

1.3 Overview of Phase II Project Objectives . 

The primary objective in Phase II was the demonstration 
of the integration of remote sensing technology with existing 
techniques for producing a vegetation type map and vegetation 
productivity estimates. A parallel objective was the transfer 
of this technology to BLM personnel to promote the implementa- 
tion and utilization of the procedures and techniques within 
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the BLM operations framevfork. The technology being deiaonstrated 
emphasized the Integration of quantitatively-based remote 
sensing data while maintaining consideration of cost-efficiency 
in the implementation of the overall design. The technology 
transfer expanded on this emphasis with particular attention 
on further developing the understanding, on the part of BLN, 
of approaches to inventories that Integrate multiple data 
sources given various resource information objectives. Specif- 
ically, any considerations based on the candidate approaches 
were focused on sampling strategies and analyses of the costs 
of data collection as specified by those strategies. 

1.4 Overview of the Approach . 

1.4.1 Technology Demonstration . 

The technology demonstration aspect of Phase II 
combined data from Landsat, low-altitude color aerial photography 
and ground visits to produce the vegetation type map as well 
as the productivity estimates for range, woodland and forest 
resources. One Landsat scene covering the Arizona Test Site 
was selected for computer-aided processing to provide the 
basis of the type map and to control the productivity sampling 
efforts (Section 2.1). Photo- interpretation results from large 
scale (1:750 - 1:1200 scale) color aerial photography was used 
to describe the vegetation characteristics of the spectral 
classes derived from the Landsat processing (Section 2.2). Also, 
estimates of parameters related to productivity were made from 
the aerial photography to be combined in a double-sampling 
freunework with ground measurements of productivity. Photo data - 
and ground data were combined through specially designed formulae 
to produce the necessary productivity estimates and also the 
vegetation species composition descriptions for each spectral 
class (Section 2.3). Map-type outputs were created based on 
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the Landsat data and the related class descriptions and, for 
some maps, productivity information was also included 
(Section 2.4). 

1.4.2 Technology Transfer . 

The technology transfer portion of Phase II was 
divided into three parts consisting of a planning session, 
workshops (two) and project status reviews (four) . The planning 
session was held at the outset of the project to familiarize 
project participants with the full scope of the work to be 
performed and to establish a milestone schedule (Section 3.1). 
The two workshops were held to provide "hands-on” instruction 
for project participants in specific Phase II technology areas 
and utilized actual project data (Section 3.2). The four 
project status reviews were held periodically through the life 
of the project to maintain continuity in reporting on interim 
results, to resolve any problems that occurred and to update 
the project schedule as necessary (Section 3.3). 

2.0 TECHNOLOGY DEMONSTRATION. 


2.1 Landsat Data Processing. 

2.1.1 Objectives . 

Given the technology objective of a vegetation type map 
and productivity estimates for range, woodland and forest resources, 
computer-aided digital classification of Landsat satellite data 
was selected as a major source of input data. Specifically, the 
classification results would provide the basis for the vegetation 
type map as well as the sampling frame for the productivity esti- 
mation tasks. The sample allocation, saiq>le selections and 
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tabulations of ^istlmate' summarlzationa all would usa tha Landsat 
classification. 

2.1.2 Procadures . 

Tha Landsat data processing for this projact was dividad 
into six major categories: 


1) 

Preprocessing 

2) 

Training 

3) 

Cla'U'if ication 

4) 

Stratification 

5) 

wiar /sscription 

i) 

Digital Terrain Data 


Preprocessing . Landsat digital data requires a certain 
anssunt it manipulation before it can be used effectively and 
effi'^iently in a resoxirca assessment/lnvantory projact. To be 
used more effectively, the data is placed into tha proper format# 
corrected for radiometric and geometric distortions# and related 
to a preselected groxind coordinate system. To be used more 
efficiently# the project area is extracted from the appropriate 
Landsat scenes to eliminate unnecessary manipulation of data that 
is outside project boundaries. 

Training . Unsupervised and supervised training 
procedures were utilized to generate the spectral statistics to 
be applied in the maximum likelihood classification algorithm used 
to classify the Landsat data. The unsupervised approach was 
used first to develop an initial set of training statistics. 
Preliminary classification based on those statistics performed 
on four subsections of the project area was evaluated for complete- 
ness. Based on that evaluation# specific training sites were 
selected and added to the original set of spectral statistics in 


2.1.2 


— Contlnu«4 


order to Improve the results from the prellminery elassiflcetion. 
The overall training procedures generated 83 spectral clusters 
(distinctly separable spectral responses identified in the 
Landsat data) to be used in classifying the Landsat data. 2t 
is important to nj>te that at this point there was only a general 
vegetation type description associated with each cluster. The 
subsequent sampllng/photo-interpretation procedures (Section 2.2) 
developed detailed quantitatively based vegetation descriptions. 

Classification . The 83 spectral clusters were applied 
in a maximum likelihood classification algorithm to the raw 
Landsat data covering the project area. The output of this step 
was a computer classification with 83 spectral classes generally 
related to the vegetation of the area. These classes were evalu- 
ated by BLM-Arizona personnel for potential of any spectral class 
representing more than one vegetation type. After all such 
"confusions” were noted, it was determined that an environmental 
stratification could be used to eliminate them. 

Stratification . Two types of stratification were applied 
to the Landsat data: administrative and environmental. The 

administrative stratification partitioned the project area into 
selected pastures and grazing allotments. This would be used to 
control the sampling and data sv'.u>;arization of the productivity 
estimation tasks. The environmental stratification divided the 
area into low desert, high desert, mountain, and forest types. 

This was used to assign new class identifiers to those spectral 
classes occurring in more than one environmental strata. 

Class Description . In order to sample the classification 
most efficiently and to reduce the total number of sa8g>les re- 
quired, it was necessary to progress from the general vegetation 
descriptions of the detailed spectral classes to preliminary 
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v«g«tation ccKnposition •■tlnatas for a aat of aunmary olaaaaa. 

Basad on tha ganaral vagatotlon typa notad for aach datailad 
■paotral data, apactral aimilaritlaa batwaan olaaaaa and proxiadty 
of locationa batwaan olaaaaa within tha anvlronnantal atrata; tha 
datailad apaotral olaaaaa wara groupad Into a aat of aunmary 
olaaaaa for aubaaquant aampllng and aatlmatlon. Thaaa aunmary 
olaaaaa wara than comblnad with axlating photo intarpratatlon 
data ovar tha projact araa in a ona-way analyala of varianca to 
produca aatimataa of tha avaraga parcant oompoalticn of aach 
claaa In tarma of traaa, ahruba, graaa, and non-vagatatlon. Tha 
summary class dascriptlons wara major inputs In datarmlnlng sampla 
slza, allocation and salaction. 

Digital Tarraln Data . Tha final a tap in tha Landsat 
procassing phasa was to input and ragiatar digital tarraln data 
with tha Z4mdsat data for usa in tha vagatation mapping task. Tha 
datailad dascriptlons of apactral classas to ba ganaratad follow- 
ing tha sang>ling of tha abova sunraary elaasas %#ould ba augsMntad 
by topographic information. Spacifically, for aach claaa tha 
rangas of alavation, slopa and aspact would ba datarminad and 
than includad with tha spacias ccxnpoaitlon dascriptlons. This 
tarraln information is raportad so to furthar axplain any apactral 
class diffarances dua to topographic Influancas. 

2.1.3 Landsat Procassing Rasults . 

Tha training staps produced 83 spectral clusters 
(training statistics) to ba ucad in classifying tha raw Landsat 
data. Elimination of "confusion" classas utilising tha anvlron- 
mantal stratification resulted in expanding to 117 datailad 
classas fr<m tha original 83 classas. Class description staps 
grouped tha 117 classes into 27 summary catagorias as shown in 
Table 2-1. Tha sections following describe tha sajqpling« 


Tebltt 2-1. lUtulti of Assignaiont of All Spoctral 
CluBtors to Summary Catagorlaa. 


I Sunimary 




Strata 


Daacrlption 


Spactral CXuatara 


Low Dasart 

Craosota-Bursaga 
(rocky soil) 

l,2,7,10ai»87,88, 

89,90,91,92,93,98, 

100 

Low Dasart 

Craosota-Bursaga 
(sandy soil) 

3,5,9,86,94,95,99 

Low Dasart 

Craosota-Pura 

4 

Low Dasart 

i 

i 

Upland Dasart 
Shrub - Craosota 
Dominant 

6,12,13,14,96,97 

Low Dasart 

Blackbrush 

45 

Low Dasart 

i 

1 

Mixad Dasart 
Shrub - Craosota 
and Cactus 

8,74 

All 

Riparian 

Woodland 

25,26,27 

High Dasart 

Shrub-Grass 

28,34,38 

High Dasart 

Grassland-Shrub 

31,39 

High Dasart 

Snakawaad-Grass 

29,30,106,107,108 

High Dasart 

Saga-Mix Shrub 

32,44,102,103,105, 

111 

High Dasart 

1 

Saga 

35,49,51,54,55, 
110,112 1 

High Dasart 

Saltshrub 

36,37,104,109 

High Dasart 

Pinyon-Junipar- 

Saga 

33,50 

High Dasart 

Pinyon-Junipar- 

Shrub 

53,56,57,58,59,61, 

62,63,84,85 

High Dasart 

Pinyon-Junipar 

40,41,46,52 

Mountain 

Mountain Shrub 

42 

Mountain 

Mixad Chaparral 

60,113 

Low Dasart 

Agricultura 

15,16,17,18,19,20, 
21,22,23,24,114 | 

Mountain 

Pondarosa Pina- 
Oak 

66 
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Strata 

Dwteription 

Spactral Cluttara 

21 

Mountain 

Pondarota Plna- 
Mia 

65 

22 

Mountain 

Pondarota Pina 

64 

23 

All 

Shadow 

47,83 

24 

All 

Watar 

67,68,69,70,71,72, 

73,82 

25 

All 

Bara 

76,77,78,79,80,81 

26 

Low Oasart 

Upland Daaart 
Shrub-Blackbruth 

43,48,75 

27 

Hijh Da tart 

Saga-Grata 

101,115,116,117 









2.1.3 


--Continued 


estimation and vegetation mapping phases and their utilization 
of the digitized administrative boundaries and digital terrain 
data. 


2.2 Data Collection. 


2.2.1 Overview of Sample Design. 

The general seunple design used in this project could 
be described as t^«o-stage with double sampling. The first stage 
was a rectangular grid of primary sanple units (PSUs) with 
dimensions of 2000 meters (east-west) by 450 meters (north-south) . 

A set of those PSUs was selected for coverage with large scale 
(1:750 to 1:1200 nominal scale) aerial photography (LSP) from 
which to derive the PSU estimates. The second stage was a double 
sample of photo plots and ground plots from within designated PSUs. 
Rangeland, woodland and forest types were sampled independently 
of one another even though the PSUs for each were defined by a 
common grid. Other features of the design varied with the vege- 
tation type being sampled as described in the following paragraphs. 

2. 2. 1.1 Rangeland. 

The quantities to be estimated were pounds per acre 
and kilograms per hectare of palatable forage for cattle within 
ten specified allotments. Figure 2-1 illustrates the locations 
of those allotments within the project area. The BLM-specif ied 
precision requirement for the estimates was as follows: the 

allowable sampling error for the overall forage estimates was 
±20% at the 80% confidence level. The results were to be tabu- 
lated by vegetation strata (as defined by Landsat processing) and 
also by pasture and allotment. The primary results were to be 
total available palatable forage for cattle adjusted for utili- 
zation. The secondary results were to be current palatable forage 
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1. Littlefield Community 

2. Blackrock 

3. Whiterock/Soapstone 

4. Wolf hole Canyon 

5. Lower Hurricane 

6 . Toquer Tank 

7. Mainstreet 

8. Jump Canyon 

9. Blackwillow-Tassi Spring 

10. Parashaunt 


Figure 2-1. Range Allotment Boundaries. 
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wh€th«r available or not. The vegetation atratun eatimatea were 
to have 80% confidence intervala. For the total eatimatea acroaa 
all atrata* the preciaion requirement would be met if one-half 
the length of the 80% confidence Interval waa no more than 20% 
of the actual eatimate. 

To allocate and aelect aamplea, the aample deaign went 
through an optimisation proceaa. A Survey Planning Model (SPM) 
developed at the University of California by the Remote Sensing 
Research Program applied summary class description data (Section 
2.1.3 Landsat Processing Results) to create a model representative 
of the magnitude, variability, and spatial distribution of forage 
from grasses, shrubs and forbs for 14 vegetation strata (1, 2, 

3, 6, 8, 17/18, 9/12/27, 10/13, 11, 14, 15, 16, 21, 26 from 
Table 2-1) . The SPM combined this data with cost information 
for collection and analysis of the photo and ground data as 
well as the 20%/80% precision requirement by means of a non- 
linear programming procedure to optimize the sampling method, 
number of samples and allocation to strata. The result was a 
stratified two-stage design with Landsat providing the strati- 
fication. The first stage was to select 108 primary sample 
units of size 450 meters by 2000 meters to be used as locations 
of the second stage sampling units. These PSUs were selected 
with equal probability from within the 14 vegetation strata. 

The second stage was in itself a double sample of aerial photo 
plots and ground plots. Each of the 108 PSUs was to be flown 
with large scale aerial photography such that each was composed 
of 15 photo plots at equal intervals. This resulted in a total 
of 1620 photo plots which were to then be photo interpreted for 
assignment into cover type strata (1447 fell within the strata 
of interest) and 135 sample plots were selected from the 14 
vegetation strata of interest for ground measurement. Each 
ground plot was to be subsampled by five parallel transects of 
40 subplots each for a total of 200 points per ground plot. 
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2. 2. 1.2 Woodland. 


Th« quantities of interest were cubic foot volume per 
acre and cubic meter volume per hectare for pinyon pine and 
juniper. Sise (diameter) and age were also of interest. The 
precision requirement was the same as for rangeland: estimates 

within ±20% at the 80% confidence level. The results were to be' 
tabulated by allotment (see Figure 2-2) and 80% confidence 
intervals were to be supplied for each. For the total estimates 
across all three allotments, the precision requirement would be 
met if one-half the length of the 80% confidence interval was no 
more than 20% of the actual estimate. 

A sample design optimization process was used similar 
to the rangeland case to specify sampling method, number of samples 
and allocation to strata. The SPM created a model representative 
of the magnitude, variability and spatial distribution of pinyon/ 
juniper volume for the three woodland strata identified and 
described during the Landsat processing phase. Based on the 
volume model, the costs associated with data collection and the 
precision requirements, a two-stage design was again specified. 

The first stage was to select 45 primary sample units with 
probability of selection proportional to the acreage of pinyon/ 
juniper within each PSU. In contrast, rangeland PSUs vrere 
selected with equal probability. The second stage was again a 
double sample of photo plots and ground plots. With 15 equally 
spaced photo plots per flight line and one flight line per PSU, 
a total of 675 photo plots were to be available from which to 
select ground plots. Each photo plot was assigned to a cover 
type strata and samples for ground measurement were selected with 
equal probabilities from those plots falling in the strata of 
interest: 





Figure 2-2. Woodland Allotment Boundaries. 
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Photo Stratiu n 

Plnyon/Junlper/Sage 
Pinyon/ Juniper/Shrub 
Plnyon/ Juniper 
Total 


No. of Photo Plots 

232 

104 

155 

491* 


No. of Ground Plot a 

5 

6 
5 

16 

the total 675 photo 


*Reflects that the remaining 184 plots of 
plots did not fall in the above strata. 


There was also to be subsampling within the photo and 
ground plots using the line intersect method to select trees to 
be measured for volume estimation purposes. One of several tem- 
plates (photo scale dependent) each with two parallel lines 
representing 75 feet each was to be overlayed on the photo plot. 

Any tree whose crown intersected a line would be photo interpreted 
for volume. For the 16 plots to be ground visited as well, the 
same photo interpreted trees were measured for volume on the ground. 


2. 2. 1.3 Forest. 

The quanitities of interest were ponderosa pine board 
foot volume per acre by stand (vegetation strata) for each area 
illustrated in Figure 2-3. The precision requirement, as in the 
previous cases, was ±.20% at the 80% confidence level. The require- 
ment would be met if for the combined areas one-half the length 
of the 80% confidence interval was no more than 20% of the actual 
estimate. 


The sampling method, number of samples and their allo- 
cation were derived by the sample optimization process described 
for the woodland case above but used summary data for three pon- 
derosa pine strata. The result was a stratified two-stage design 
where the first stage was to select 47 PSUs with equal probabili- 
ty from within the three strata. (In contrast, the %#oodland PSUs 
were selected with probability proportional to area.) 
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The second stage was a double sample of aerial photo 
plots and ground plots. Each PSU was covered with 15 photo plots 
at equal intervals for a total of 690 photo plots (one PSU photo 
set was unusable resulting in 15 less plots than anticipated) . ' 

The photo plots were interpreted to be assigned to cover type 
strata and semples were selected from the ponderosa pine strata 
for ground measurements: 

t of PSUs I of PSUs I of Photo Plots f of Ground 
Stratum Selected Available within strata Plots to Visit 

Ponderosa Pine/ 


Oak 

4 

4 

56 

4 

Ponderosa Pine/ 
Mixed Shrub 

26 

25 

62 

15 

Ponderosa Pine 

17 

17 

47 

10 

Total 

47 

46 

165* 

29 

*Reflects that 
ponderosa pine 

525 Of the 
strata. 

690 photo plots 

were not assigned 

to 


There was subsampling within the phot < and ground plots 
using the line intersect method as described for the woodland case 
in Section 2. 2. 1.2 above. 


2.2.2 Summary of Allocation and Selection of Samples. 

Results of the procedures described in the sample design 


overview above are summarized in the 

t of Strata I of PSUs 
Vegetation of Interest to Select* 

following table: 

V of Photo Plots 
Assigned to 
Strata of Interest 

1 of 

Ground Plots 
Selected 

Rangeland 

14 

108 

1447 

136 

Woodland 

3 

45 

491 

16 

Forest 


46 

165 

M 

Total 



20 

199 

2103 

181 

*Assumes loss 

of one PSU 

in forest type. 
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2.2.2 —Continued . 

The number of photo plots above reflects that out of 
2985 possible plots, only 2103 were assigned to strata to be 
sampled for productivity. All photo plots were interpreted for 
vegetation composition in producing the vegetation map of the 
area. Of the 181 ground plots selected, two forest plots did not 
contain any trees and therefore were not visited on the ground; 
and one range plot was on such steep ground that measurement was 
Impossible. As a result, 178 ground plots were actually visited 
and the three missing plots were not replaced as this could have 
introduced unwanted bias into the productivity estimates. 

2.2.3 Photointerpretation on the Large-Scale Aerial Photography . 

As described above, each photo plot from the PSUs was 
initially interpreted to assign it to a cover type strata as part 
of the procedures to select second stage samples. There were, 
however, two more phases of photo interpretation: vegetation 

composition and productivity estimation. Each photo plot was 
first interepreted, independently of cover type assignment, for 
percent species composition and then used to describe the Landsat 
spectral classes. Finally, photo plots within the forest and wood- 
land vegetation types were interpreted for tree volxime estimates 
then used to estimate volumes for the Landsat summary classes 
(cover-type strata) . These procedures are described in more 
detail below. 

2. 2. 3.1 Vegetation Composition. 

The large scale aerial photography (LSP) consisted of 
35mm format color slides flown to achieve nominal scales ranging 
from 1:750 to 1:1200 (actual scales ranged from 1:500 to 1:3000 
due to rugged terrain in some areas) . An LSP photo plot was a 
stereo pair with the stereo overlap defined as the area of in- 
terest. There were 15 equally spaced plots along the center of 
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the long dimension of each PSU (2000 meters by 450 meters) 
resulting in plot centers approximately 125 meters apart and each 
plot representing 60 feet to 100 feet square on the ground, 
depending on photo scale. 

To combine the photo plot Interpretation with the 
Landsat spectral classes, the location of each plot was Identified 
on uses quads and then digitized and transformed to the same 
ground coordinate system to which the Landsat data had been 
registered. This allowed matching of the Landsat classes to the 
interpreted data on a polnt-by-polnt basis. 

Within the stereo overlap of each photo plot the follow- 
ing vegetation composition Information was interpreted and 
recorded: identification of major plant species occurring on 

the plot, percent ground cover for each species, and also height 
class and foliar density by species. The Interpretation «ras per- 
formed by one person to eliminate multiple biases that would occur 
by using more than one Interpreter. Thus, any bias In the Inter- 
pretations of the plot characteristics was constant. 

2. 2. 3. 2 Woodland and Forest Volume. 

In addition to vegetation composition, trees were 
selected via the line Intersect method to interpret for volume 
estimation on all plots from PSUs covering the forest and wood- 
land vegetation types ( 676 photo plots out of the total 2985) . 

One of several templates (depending on the plot photo scale) was 
overlayed in a predetermined manner (to ensure randomness) on to 
the stereo overlap portion of each photo plot. Each template had 
two parallel lines with each line representing 75 feet on the 
ground at the photo scale for which the ten^late was prepared (a 
total of nine templates were prepared representing photo scales 
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2. 2. 3. 2 —Continued 


from It 500 to 1x3000 to accommodate the range of actual photo 
scales noted in the LSP) . For a given plot, any tree whose crown 
touched one of the lines, information wcs determined and recorded 
as to its species (pinyon pine, juniper o;r ponderosa pine were 
the only species of interest) , crown diameter (measured) and 
height (estimated to the nearest five feet) . 

2.2.4 Ground Data Collection. 

The data collection procedures for the ground plots 
were designed to be as close to BLM procedures as the multistage 
design would allow. The procedures used on the range plots were 
modified from the Soil Vegetation Inventory Method (SVIM) while 
the forest and woodland procedures were based on the BLM Extensive 
Forest Inventory Field Handbook. 

2. 2. 4.1 Range Plots. 

Plot location was defined as the center point of the 
stereo overlap of a selected LSP plot. After locating it on the 
ground, this point was used as the center of a 200 point rectangular 
grid as formed by 5 lines with 40 points per line spaced evenly 
over a 42 foot by 48.75 foot area. These 200 points were used in 
determining species composition of the plot. A systematic set 
from the 200 points of 20 subplots (.1 square meter) were used to 
take SVIM shrub characterization and ocular weight estimates, 
recorded on modified SVIM field data forms. A subset of the ocular 
weight estimate plots were clipped, air-dried and weighed to use 
for adjusting the ocular estimates. This was done for each field 
crew member so as to have Individually developed adjustments for a 
given plot based on who made the original ocular estimates. 
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2, 2. 4, 2 Woodland and ror««t Plots. 

Ai described previously, plot location was determined 
by the parallel transects in the stereo overlap of each selected 
LSP plot. These transects %irere annotated on the LSP photos to 
be used in the field with the individual trees to be measured, 
marked and numbered to ensure correspondence with the photo 
interpretation data. Por each tree, the following information 
was recorded: diameter (ground and stump height for woodland and 

breast height for forest) , height, crown diameter, number of stems 
per tree (woodland) and age (on a subset of trees) . Measurement 
of the distance between two photo>ldentified ground-located 
points was also recorded for use later in determining actual 
photo scale. 

2.3 Productivity Estimation and Map Verification. 

The following sections describe the steps and procedures 
used to combine tne ground measurements, photo estimates and 
Landsat classification to produce the outputs specified by the 
project objectives. 

2.3.1 Vegetation Type Mapping, 

Before combining the Landsat classification results with 
the photo interpretation data (species composition) , the 2985 
photo plots were separated into two groups. The first group was 
composed of 2399 plots which were used to generate the vegetation 
descriptions of the Landsat classes. The remaining 586 photo 
plots, the second group, were used to verify those descriptions. 
This partitioning of photo plots was done in such a manner as to 
ensure that an adequate number of photo plots would be available 
to both describe and then verify the maximum number of Landsat 
spectral classes. 
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2.3.1 — Contlnu»d . 

Th« photo data from each plot from tha firat group was 
than natchad with tha Landsat classification on a point-by-point 
class-by-clasa basis using analysis of variance tachniquas. 

This produced tha mean and standard deviation of percent cover 
by species for each Landsat class. The same procedures were 
followed with the verification set of plots. BLM-Arisona personnel* 
using tha vegetation descriptions for the 117 Landsat classes and 
checking locations of the classes within the project area* as- 
signed each class to a category within the Arizona vegetation 
framework. This resulted in 14 Level III vegetation framework 
categories which could be aggregated further into 10 Level II 
categories. Table 2-2 shows the results of the assignments. The 
digital terrain information, elevation* slope and aspect was com- 
bined with the framevrork categories resulting from above to augment 
the species composition descriptions of each. Appendix II con- 
tains the "menu” descriptions (vegetation and topography) of each 
of the 14 Level III categories. The verification data set was 
compared with the original description set and it was noted that 
no significant differences existed between them. The descriptions 
of the vegetation framework classes and their locations (based on 
the Landsat data) were therefore considered accurate representations 
of the ground conditions in the project area. 

Area estimates by vegetation type with 80 percent 
confidence intervals at the pasture and allotment level were 
generated by relating the Landsat class counts to their descrip- 
tions from the photo data through linear regression. In the 
example in Table 2-3* the ”ad'ju»ted estimate" represents the 
regression coefficients applieu to the Landsat area summaries. 

The ”80% confidence" column represents an Interval which has ap- 
proximately 80% probability of containing the actual value. 
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Tftbl* 2-2. Aisignm«nt of th« 117 Computor Clatsos 
to th« Vagatation Pramawork. 


Vagatatlon Pratnawork Claaaif icatlon Cowputar Claaaaa - 


1) 

1— 

Agrieultura 

15-24 « 114 

2) 

211 

Pondaroaa Pina Porast 

64, 65 

3) 

311 

Pinyon-Junipar Koodland 

33, 35, 41-43, 46, 47, 
52, 53, 55-59, 61-63, 
66, I3-B5 

4) 

322 

Riparian Hoodland 

2j, 26 

5) 

412 

Upland Daaart Shrub 

1-10, 12-14, 74, 75, 
81, 86, 87, 89-100 

6) 

421 

Graat Basin Sagabrush 

11, 28-30, 32, 34, 38, 
40, 44, 48-50, 54, 
101-103, 105, 106, 108 
no, 112, 115-117 

7) 

423 

Blackbrush 

88 

8) 

424 

Othar Tall Shrub 

111 

9) 

425 

Half Shrub 

45 

10) 

432 

Oakbrush 

60, 113 

11) 

433 

Othar Mountain Shrub 

27, 31, 39, 51 

12) 

511 

Parannial Grassland 

36, 37, 77, 10a, 107, 
109 

13) 

6— 

Barran Land • 

76, 78-80 

14) 

7— 

Watar 

67-73, 82 
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Table 2-3. Example of Area Estimates in Hectares 
by Allotment and Pasture. 


. Vegetation Type: Pinyon/ Juniper Woodland 


ALLOTMENT 

PASTURE 

LANDSAT-ONLY 

ESTIMATE 

ADJUSTED 

ESTIMATE 

80% CONFIDENCE 

LOWER HURRICANE 

489 

494 

285,703 ' 

ST. GEORGE 

192 

* 

194 

112,276 

COYOTE 

109 

no 

63,156 

NORTH GYP 

75 

76 

44,108 

WEANING 

3 

3 

2,4 

SOUTH GYP 

23 

23 

13,33 

GRAVEL 

60 

60 

35,86 

PETE 

23 

23 

13,33 

HOLDING PEN 

0 

0 

0,0 

PARKER 

1 

1 

1,1 

FORCE 

5 

5 

3,6 

SELLING 

0 

0 

0,0 

TOQUER TANK 

44 

44 

26,63 

1 

5 

5 

3,7 

2 

39 

40 

23,56 

MAINSTREET 

6807 

6875 

3971,9779 

CECIL 

7 

6 

4,9 

ROUKT POND 

1 

1 

0,1 

SQUARE POND 

8 

8 

5,11 

ANTHONY'S HIGLEY 

16 

16 

9,23 

CALVING 

12 

12 

7,17 

TEMPLE TRAIL 

27 

27 

16,39 

SALARATUS 

63 

64 

37,91 

MUDHOLE 

19 

19 

11,27 

TWIN TANKS 

367 

371 

214,527. 

WARDS 

23 

23 

14,33 

DUTCHMAN 

693 

699 

404,995 

COX-ATKXN 

1511 

1526 

881,2171 

ENGLESTEAD 

1973 

1992 

1151,2834 

LITTLE JOE 

2036 

2057 

1188,2925 

BISHOP t BURR 

53 

53 

31,76 
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2.3.2 


Rangeland Productivity Estinvation. 


The data reduction process started with the air-dry 
weights of the species clipped on the various subplots of each 
range ground plot (see Section 2. 2. 4.1) and related these to the 
corresponding ocular estimates through linear regression (R^“77%), 
The resulting regression coefficients were then used to adjust 
the remaining ocular subplot weight estimates to produce current 
forage weight estimates, "available" or not, the secondary 
estimates. At t. . > point, the availability and utilization 
factors were applied to produce the primary estimates; forage 
weights "available" whether "utilized" or not. Both types of 
estimates were carried through the remainder of the estimation 
process; however, at this stage all estimates for unpalatable 
species (for cattle) were deleted from the process. The adjusted 
estimates were averaged by species, weighted (statistically 
speaking) by the frequency of species occurrence in all subplots, 
and then summed to obtain forage weight estimates acrosc« all 
palatable species for each of the 135 ground plots. 

As described previously (Section 2. 2. 1.1), the LSP was 
used as a stratification tool in selecting the plots to be 
visited on the ground. Attempts were made to correlate the 
ground-based estimates of forage weight with photo-interpreted 
values such as crown cover, the cover/density product, and the 
cover /density/height product, but the relationships were not 
strong enough to be useful. Consequently, for the range case, 
the LSP was used only for stratification and not estimation. 

The ground-based weight estimates were related to their 
corresponding Landsat summary classes to produce weight estimates 
per unit area by class, including the associated 80% confidence 
inter yal for each. Forage weight estimates for the specified 
allotments were developed by ratio estimation. (See Appendix III 
for summary of estimates.) 
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2.3.3 


Woodland Volume Estimation. 


The ground-based measurements of diameter at stump 
height, total tree height and average crown diameter for trees 
selected using the line Intersect method were converted to cubic 
foot volumes through BLM pinyon- juniper volume tables. The 
ground volumes were then related to the Photo values through 
multiple regression (r 2 ■ 37% for pinyon? 34% for juniper) to 
yield a tree volume prediction equation which could be used for 
the photo-interpreted trees. These predictions were then 
expanded through an estimation procedure that used the tree crown 
diameters (which were proportional to the probability of selection 
for a tree) and the transect lengths to develop volume per acre 
estimates for each photo plot. These were then averaged by PSU 
and combined with PSU sampling weights to obtain the final cubic 
foot volume estimates (also converted to metric measure: cubic 
meters per hectare) for the three allotments of Interest. 

Average diameter and age estimates were obtained from ground 
measurements alone. Estimates of precision were also produced 
using methods fitting the sample design and parameter estimation 
process. (See Appendix IV for summary of estimates.) 

2.3.4 Fo rest Volume Estimation. 

The ground-based measurements of diameter breast height 
and total tree height for trees selected by the line intersect 
method were measured and used to compuv board foot volume using 
equations generated by the U.S.* Forest Service (BLM-Arirona 
personnel designated these equations as acceptable ones to use 
for the area) . These ground determined volumes were then re- 
gressed against the photo-interpreted values (height and crown 
diameter) for the same trees (r 2 ■ 78% for ponderosa pine) to yield 
a tree volume prediction equation which could be used for all the 
photo-interpreted trees. Both the photo and ground tree volume 
estimates were expanded to board foot per acre estimates for the 
entire plot by use of the crown diameters (which were proportional 
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2.3.4 


— Continued. 


to the probability of selection for a tree) and transect lengths. 
The matched photo and ground plot estimates were then related 
through regression for the Landsat strata sampled (R^ ■ 44% for 
ponderosa pine/oak and ponder osa pine/mixed shrub; • 76% for 
ponderosa pine pure) . A statifled regression estimator was then- 
used to obtain the final estimates by Landsat strata, and the 
associated precision estimates were also produced. (See 
Appendix V for summary of estimates.) 

2.4 Costs Analysis. 

The costs associated with this project were separated 
into recurring and non-recurring components. Only recurring 
costs were considered in this analysis. In other words, "false" 
starts, development of procedures, and non-special management 
costs were omitted to the extent possible. The objective of this 
analysis was to associate a cost per acre with both the vegeta- 
tion mapping and the productivity estimation elements of this 
project. Table 2-4 summarizes rhis analysis and shows that at 
the average price per ESL manhour over the life of the project, 
the mapping effort for the 2.2 million acre site cost 3^ per acre. 
The productivity estimation, performed on 550,000 acres, cost 
16.3^ per acre given the mapping had been performed. These figures 
were based on the costs associated with achieving the objectives 
and accuracies of this project. If similar projects were under- 
taken but with higher accuracy and precision specifications, it 
should be anticipated that the costs per acre could be greater. 

2 . 5 Landsat Digital Products. 

The Landsat digital hardcopy output products consisted 
of 4 application maps and 2 sets of color-coded classification 
images. The application maps were generated from the digital 
classification quantitative inventory and ancillary data in the 
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Table 2-4. Phase II Cost Summary. 


Project Element 

Total ESL , 
Man-Hours 

Mapping 

Component 

^ Estimation 
Component 

Additional 

Costs 

Landsat Processing 

570 

450 

120 


Multistage Sampling 

270 

120 

150 

4800^ 

Data Collection 
Photo Acquisition 

50 

50 


11,300^ 

Photo Interpretation 

700 

500 

200 


Ground Data Collection 

100 

mi 

100 

31,970^ 

Productivity Estimation 

1250 

240 

1010 

700^ 

4 

Output Products 

Map Output 

180 

180 

- 

, 

Tabular Summaries 

120 

40 

80 


Totals 

3240 

1580 

1660 

$48,770 


Vegetation mapping cost ■ [(1580 hrs x $3. 17/hr.) + 11,700] 

* 2.2 million acres ■ 3(/acre 

Productivity estimation « (Given the mapping is completed and avail- 
able) 

- [(1660 hrs X $31. 70/hr) + 37,470j 
i 550,000 acres « 16.3C/acre 


1. R5RP, U.C. Berlceley, in survey planning model runs and estimation 
worlc. 

2. Acquisition of 3000 large scale aerial photography plots. 

3. Collection of ground data on 179 ground plots. 

4. Does not include generation of final report. 
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2.5 


— Continued . 


form of digital terrain data and administrative data. The map 
products were: 

a. Rangeland Suitability, defined as: 

(1) current production of useable forage above 
20 lbs. per acre, 

(2) located within a four mile radius of water, 

(3) slope less than 51 percent. 

b. Potential Rangeland Suitability, defined as: 

(1) current production of useable forage above 
20 lbs. per acre, 

(2) service area of water is greater than 
four mile radius, 

(3) slope less than 51 percent. 

c. Sagebrush Treatment Area, defined as: 

(1) Great Basin sagebrush vegetation community, 

(2) slope less than 15 percent. 

d. Fire Flash Fuel Areas, defined as: 

(1) annual grass and forbs cover 30 percent 
or greater of ground, 

(2) aspect southwest - south - southeast. 

These products were presented with a mapping minimum of 10 acreas 
(4 hectares) at a scale of 1:126,720 as a halftone grey-coded 
map on stable base transparency material. Maps 1 and 2 were 
limited to the 10 allotments sampled for range productivity. 

Maps 3 and 4 represent data over the full project area. 
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2.5 


— Continued. 


The 5th map output products were for the forest admin- 
istrative area A and consisted of 8x10 inch Polaroid hardcopies 
of the color-coded computer class groupings. These products 
were presented at 10 acre (4 hectare) minimum mapping units for 
the Level II framework and for the Level III framework groupings 
developed by the BLM-Arizona personnel. An IDIMS computer 
compatible tape (CCT) was generated with this data for the full 
project area as well. 

3.0 TECHNOLOGY TRANSFER ACTIVITIES. 

3 . 1 Planning Session - October 3-4, 1978. 

The objective of this session was to familiarize all 
project participants with the expected flow of tasks over the 
life of the Phase II efforts. The project flow chart in Appendix 
I was presented at this time. Also, the two training courses 
were discussed and finalized as to scope, content and dates. 

A preliminary set of milestone dates were established to be used 
in monitoring progress on the numerous aspects of the Phase II 
effort. 

3 . 2 Seminars/"Hands-On" Training. 

There were two training courses given to BLM-Arizona 
personnel during Phase II; Multistage Sampling-Field Workshop 
and the Data Analysis Workshop. 

3.2.1 Multistc S ampling-Field Workshop, May 21-25, 1979. 

This week long workshop/seminar cov^ired the sample 
allocation, data collection and data analysis procedures used on 
the project. The objective of the session was; 


3.2.1 


— Continued 


a. The development of a thorough, practical under- 
standing of sample allocation procedures used on 

the project and the justification for the procedures. 

b. To introduce and develop, from a practical stand- 
point, the statistical procedures and estimators 
being used on the project. 

c. To further develop understanding of the 
optimization procedures being used with an em- 
phasis on the factors that must be considered 
when designing an inventory and mapping project. 

d. To further the understanding of the data sources 
being used on the project, including (1) the 
information extraction procedures, (2) relative 
level of accuracy and precision of each data 
source, (3) relative cost of acquisition and 
data extraction for each source, and (4) the 
risk associated with the use of data from each 
source. 

e. To prepare the attendees for the data analysis 
workshop to be held at ESL in Fall 1979. 

Participants in the course stated afterwards that the 
objectives were, in general, satisfied. The most frequent comment, 
however, was that of having a much greater appreciation for the 
complexity of resource inventory programs such as Phase II- 

3.2.2 Data Analysis Workshop, November 13-16, 1979. 

The objective of this course was to familiarize program 
participants with the procedures and techniques to be followed in 
reducing the data collected during Phase II to produce the 
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3.2.2 


—Continued . 


vegetation map and the productivity estimates. The course 
included presentations on analysis of variance, contingency 
analysis, regression analysis and ratio estimation procedures, 
all using actual data from Phase II. The most significant result 
of the course was the grouping of the 117 Landsat spectral 
classes into 14 Arizona vegetation framework Level III cate- 
gories performed by the attending BLM-Arizona personnel. 

3 . 3 Project Status Reviews. 

The objectives of these periodic reviews were to present 
progress to date on the Phase II effort and identify and resolve 
any problems noted or anticipated. The four in-progress reviews 
were held in January, June, September, and November, 1979. 

These reviews were considered extremely valuable in maintain- 
ing involvement in the project by all BLM, NASA and contractor 
participants. Coupled with the training course, the reviews 
provided an excellent means of exchanging information on the 
project; its objectives, the procedures, the problems encoun- 
tered and their solutions and the results. 

A final project review was held on 29-30 May 1980 to 
present the accomplishmencs and results of Phase II to a larger 
audience than that addressed at the status reviews. In the words 
of the NASA Project Manager for this APT, "Overall, I considered 
the Phase II Final Review as a successful culmination of a very 
rewarding and productive project. The success can be directly 
attributed to the cooperation and determination of the people 
involved from both BLM and the contractor." 
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APPENDIX I 


PROJECT FLOW DIAGRAM 


The following flow diagram illustrates the order of 
the detailed tasks and their Interrelationships with other tasks 
that were necessary to completely define the approach to the 
Phase II technology demonstration. 
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APPENDIX II 


COMPUTER MENU DESCRIPTIONS 
LEVEL III FRAMEWORK CLASSIFICATION 


This appendix contains selected examples of the 
computer class menu descriptions. See Volume II, 
Appendix 2-G for the complete set. 
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original page 18 

OF POOR QUALITY 

NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Sunmary Class: 3 No. of Photo Sanplas: 1.222 X of Total: 58 

Name: Pinvon- Juniper No. of Acres: 1,056,976 No. of Hectares: 

Wood land 

Spectral Classes: 33,35,41-43,46-47, X of Area: -41 
52-53 , 55-59 , 61-63 , 66 , 83-85 

I. VEGETATION 
(X Cover by Species) 


A. 

Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 


Ponderosa Pine 

1.0 

.1 

Creosote 




Piny on Pine 

6.5 

.J 

Bursage 




Juniper 

17.8 

.4 

Blackbrush 

.1 

.0 


Other Tree 

1.5 

.1 

Big Sagebrush 

7.0 

.3 





Other Shrub 

2.1 

.1 

C. 

Shrubs — ^Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 


Gambel's Oak 

3.2 

.3 

Cottonwood 




Turbine 11a Oak 

1.6 

.J 

Willow 





Other Shrub 

7.0 

.4 

Other Shrub 



E. 

Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 


Perennials 

.3 

.1 

Yucca 

.1 

0 


Annuals 

2.7 

.3 

Other Cactus 



G. 

Non-Vegetation 

Mean 

Std.Err. 





Barren (Rocky) 

19.7 

.7 





B'rren (Sandy) 

29.1 

.8 





Water 








Shadow 

A 

.1 
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COMPUTER CLASS 3 (Pinyon Juniper Woodland) GEOBLOCK Arizona 



ASPECT (FLAT) 






ORIGINAL PAGE IS 
OF POOR QUALITY 

NASA/BLM APT PHASE II 

ARIZONA CLASSIFICATION VEGETATION/TERRAIN MENU DESCRIPTION 

Suoanary 6 No. of Photo Samples: 535 X of Total: 25 

Name: Great Basin Sage- No. of Acres: 904,389 No. of Hectares: 361,756 

Spectral Classes 3. 1, 28-30, 32, 34. 38, X of Area: 

40, 44,48-50, 54,101-103,105-106, 108, 110,112,115-117 

I. VEGETATION 
(X Cover by Species) 


A. Trees 

Mean 

Std.Err. 

B. Shrubs — Desert 

Mean 

Std.Err 

Pondcrosa Pine 



Creosote 

.1 

0 

Plnyon Pine 

.2 

0 

Bursage 

.3 

.1 

Juniper 

1.9 

.2 

Blackbrush 

.4 

.2 

Other Tree 



Big Sagebrush 

12.3 

.7 




Other Shrub 

8.7 

.4 

C. Shrubs — Mountain 

Mean 

Std.Err. 

D. Riparian Woodland 

Mean 

Std.Err 

Gambel’s Oak 



Cottonwood 



Turbinella Oak 

.2 

.1 

Willow 



Other Shrub 

.8 

.2 

Other Shrub 



E. Grasses 

Mean 

Std.Err. 

F. Cactus 

Mean 

Std.Err 

Perenn lals 

3,2 

.5 

Yucca 

1.5 

.1 

Annuals 

14.9 

1.0 

Other Cactus 

.3 

.1 

G. Non-Vegetation 

Mean 

Std .Err. 




Barren (Rocky) 

27.2 

1.4 




Barren (Sandy) 

26.5 

1.2 




Water 






Shadow 

1.4 

.5 
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APPENDIX III 


onOR QUAUTT 


RESULTS OF RANGE PRODUCTIVITY ESTIMATION 


The primary estimates of palatable cattle forage are _ 
in Table 1 for Landsat strata and Table 2 for allotment and 
pasture. These are the amounts available regardless of utili- 
zation. The secondary estimates, the current amounts 
regardless of availablility are in Table 3 for Landsat strata 
and Table 4 for allotment and pasture. The Landsat strata 
listed are the summary categories used in sampling and esti- 
mation. Results are tabulated in both pounds per acre and 
kilograms per hectare. The paired values in parentheses are 
80% confidence intervals based on the estimated standard errors. 
To check the precision requirement for the primary results, it 
is seen that one-half the length of the confidence interval 
across all strata is (16 .7 - 10 . 6) /2 » 3.05 Ib/acre, which is 2?% 
of the estimate of 13.66 Ib/acre. For the secondary estimates, 
the result is 23%. So the precision requirement was not quite 
met, although it was fairly close. 

The classification results combined into the 27 
summary classes and the overlayed allotment and pasture 
boundaries provide a detailed cover type map illustrating the 
locations of the various range strata. 
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(U INSUFFICIENT SAMPLES Page \ of 3 



INSUFFICIENT 






















































RANGE FORAGE ESTIMATES BY LANDSAT STRATA 

PALATABLE SPECIES ONLY AMOUNTS AVAILABLE AND PROJECTEF FOR FULL UTILIZATION 



(I) INSUFFICIENT SAMPLES 











Table 2. 

RANGE FORAGE ESTIMATES BY PASTURE 




















































RANGE FORAGE ESTIMATES BY PASTURE 
































































RANGE FORAGE ESTIHATES BY PASTURE 






























RANGE FORAGE ESTIMATES BY PASTURE 



TOTAL 75,607 30,597 14.2 



RANGE FORAGE ESTIMATES BY PASTURE 
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RANGE FORAGE ESTIMATES BY PASTURE 

PALATABLE SPECIES ONLY - AMOUNTS AVAILABLE AND PROJECTED FOR FULL UTILIZATION 
















































RANGE FORAGE ESTIHATES BY PASTURE 

PALATABLE SPECIES ONLY - AMOUNTS AVAILABLE AND PROJECTED FOR FULL UTILIZATION 
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RANGE FORAGE ESTIMATES BY LANDSAT STRATA 


















































































































PALATABLE SPECIES ONLY CURRENT AMOUNTS AVAILABLE OR UNAVAILABLE 

















RANGE FORAGE ESTIMATES BY PASTURE 
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RANGE FORAGE ESTIHRTES BY PASTURE 



TOTAL 108,881 44,063 20.5 19.1 21.4 





































RANGE FORAGE ESTIMATES BY PASTURE 







RANGE FORAGE ESTIMATES BY PASTURE 
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APPENDIX IV 


RESULTS OF WOODLAND PRODUCTIVITY ESTIi'ATION 


The juniper estimates by allotment are in Table 1 
for standard measures and Table 2 for metric measures. The 
corresponding pinyon pine estimates are in Tables 3 and 4. 

In each table the "Total Area" is for the entire area of 
interest while the "Sampling Freune Area" is for BLM land only 
since other ownerships were excluded from the seunpling. The 
volumes per unit area are obtained by dividing the total 
volumes by the Sampling Frame Area. The paired values below 
the volume estimates are 80% confidence intervals. These are 
estimated from the standard error of the estimate and the 
expected deviation for 80% confidence based on a sampled 
Gaussian distribution. To check the precision requirement 
for juniper, it is seen on Table 1 that one-half the length of 
the confidence interval across all allotments is (185. 2-142. 3)/2 
= 21.4, which is only 13% of the estimate of 163.7 cubic feet 
per acre. Therefore, this precision requirement was met. 
However, for pinyon pine in Table 3, (37.0-19.0)/2 « 9.0 
which is 32% of the 28.0 estimate. This shows that the 
precision requirement was not met for pinyon pine. The column 
labeled "Relative Standard Error" contains the standard error 
divided by the estimate. The other columns in the tables are 
for average size (diameter) and average age in years, along 
with associated 80% confidence intervals. 
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Table !• 
JUNIPER ESTINRTES 


ORIGSMAL PASS Kj 
OF POOR QUALITY 
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NUMBERS IB FARCBTHCSES ARE 801 CONFIDENCE INTERVALS. 





Table 2. 

JUNIPER ESTIMATES IN METRIC UNITS 
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PINYON ESTIMATES 













Table 4 . 

PINYON ESTIMATES IN METRIC UNITS 











APPENDIX V 


RESULTS OF FOREST PRODUCTIVITY ESTIMATION 


The estimates of Ponderosa pine board foot volume 
by Landsat strata are given in Table 1 for the Blackrock-Virgin 
Mountains Area and Table 2 for the Mt. Trumbull-Whitmore Canyon 
Area. The combined totals are in Table 3. In each table the 
"Landsat Stratum" column refers to the 27 summary classes. 

"Total Area in Acres" is the entire area of interest and 
"Sampling Frame Area in Acres" is the set of all PSU's which 
contain any pixels in any of the three Ponderosa pine strata. 

The volume estimates are in the next two coluiruis. The values 
in parentheses are 80% confidence intervals. These are 
estimated from the standard error of the estimate and the 
expected deviation for an 80% confidence Interval based on a 
sampled Gaussian distribution. To check the precision requirement 
it is seen that one-half the length of the confidence interval 
across all strata is (147-113)/2«17, which is only 13% of 
the estimate of 130 board feet per acre. Therefore, the 
precision requirement was met easily. The column labeled 
"Relative Standard Error" is the standard error divided by the 
estimate. 


V-1 









2 



























{\) NO SAMPLES. INSIGNIFICANT ACREAGE - ZERO VOLUME ASSIGNED 

(2) NO SAMPLES, SIGNIFICANT ACREAGE - POOLED ESTIMATE USED FROM ALL STRATA EXCEPT 20, 21 

(3) NOT ESTIMABLE - INSUFFICIENT SAMPLES 

OTHER NUMBERS IN PARENTHESES ARE 801 CONFIDENCE INTERVALS 
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OTHER HUH8ERS IN PARENTHESES ARE BOX CONFIDENCE INTERVALS 










































(1) NO SAMPLES, INSIGNIFICANT ACREAGE - ZERO VOLUME ASSIGNED 

(2) NO SAMPLES, SIGNIFICANT ACREAGE - POOLED ESTIMATE USED FROM ALL STRATA EXCEPT 20, 21, 22 
U) NOT ESTIMABLE - INSUFFICIENT SAMPLES 

OTHER NUMBERS IN PARENTHESES ARE 801 CONFIDENCE INTERVALS 


























































